The composition and size distribution of the dust in 
the coma of comet Hale-Bopp 



M. Min, a J. W. Hovenier, a A. de Koter, a L. B. F. M. Waters, a ' b 

C. Dominik a 

^Astronomical institute Anton Pannekoek, University of Amsterdam, Kruislaan 403, 1098 

SJ Amsterdam, The Netherlands 

h lnstituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200B, 3001 

Heverlee, Belgium 



Abstract 

We discuss the composition and size distribution of the dust in the coma of comet Hale- 
Bopp. We do this using a model fit for the infrared emission measured by the infrared space 
observatory (ISO) and the measured degree of linear polarization of scattered light at var- 
ious phase angles and wavelengths. The effects of particle shape on the modeled optical 
properties of the dust grains are taken into account. Both the short wavelength (7-44 //m) 
and the long wavelength (44- 120 /mi) infrared spectrum are fitted using the same dust pa- 
rameters, as well as the degree of linear polarization at 12 different wavelengths in the 
optical to near infrared domains. We constrain our fit by forcing the abundances of the 
major rock forming chemical elements to be equal to those observed in meteorites. The 
infrared spectrum at long wavelengths reveals that large grains are needed in order to fit the 
spectral slope. The size and shape distribution we employ allows us to estimate the sizes of 
the crystalline silicates. The ratios of the strength of various forsterite features show that the 
crystalline silicate grains in Hale-Bopp must be submicron sized. On the basis of our analy- 
sis the presence of large crystalline silicate grains in the coma can be excluded. Because of 
this lack of large crystalline grains combined with the fact that we do need large amorphous 
grains to fit the emission spectrum at long wavelengths, we need only approximately 4% 
of crystalline silicates by mass (forsterite and enstatite) to reproduce the observed spectral 
features. After correcting for possible hidden crystalline material included in large amor- 
phous grains, our best estimate of the total mass fraction of crystalline material is ~7.5%, 
which is significantly lower than deduced in previous studies in which the typical derived 
crystallinity is ~20 - 30%. The implications of this low abundance of crystalline material 
on the possible origin and evolution of the comet are discussed. We conclude that the crys- 
tallinity we observe in Hale-Bopp is consistent with the production of crystalline silicates 
in the inner solar system by thermal annealing and subsequent radial mixing to the comet 
forming region (~30 AU). 
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1 Introduction 



Comet Hale-Bopp is undoubtedly the best studied long period comet in the solar 
system. The unusual brightness of this comet allowed for its discovery when it 
was still at 7 AU distance from the Earth, and a long term mo nitoring of the dust 
activity after perihelion up to almost 13 AU (TWeiler et al.Ll2003l) . The infrared spec- 
trum taken by the Infrared Space Observator y (ISO) provided a unique opportunity 
to stu dy the composition of cometary dust (|Crovisier et all 1 19971: iLellouch et all 
Il998l) . The strong resonances visible in this spectrum were attributed to the pres- 
ence of crystalline silicates. 



The relatively sharp 9.8 fim feature in the ISO spectrum together with observations 
of an extremely high degree of linear polarization lead to the conclusion that the 
dust in t he coma of Hale-Bo pp has an overabu ndance of submicron sized grains 
(see e.g. lHanneretall [T999: Mas on et all 1200 lL and references therein) 



The composition of the dust in Hale-Bopp has been modeled frequently in the liter- 
ature applying various approaches (see e.g Lellouch et all 1998 [ Li and GreenbergL 



1998 



: brucatoetaTiri 999: 



Wooden et al 



1999; 



Harker et al 



2002; 



Bouwman et al., 



2003 



Moreno et al., 2003). The interpretation of the results obtained from the 
observations lead to discussio ns in the litera ture on the origin and evolution of 
cometary dust (for a review see IWoodenL l2002f) . Cometary dust is believed to be the 
most primitive material present in the solar system. The high crystallinity found by 
most studies of the infrared emission spectrum presents a problem for this scenario. 
Since the dust that entered th e proto-solar nebula from the interstellar m edium is 
almost completely amorphous dLi and Dra ine. 2001t lKemper et all |2004). the crys- 
talline silicates in Hale-Bopp have to be products of processing in the early solar 
system. Crystalline silicates can only be produced in high temperature environ- 
ments (> 1000 K), so close to the Sun, while comets form in regions where the 
temperatures are low enough for water ice to exist (< 160 K). Thus, the crystalline 
silicates have to be transported outwards to the regions where the comets form, 
or the dust temperat ure must be loca l ly increased due to for example lightning 
or sho ck annealing (Har ker and DeschL 20021: IPilipp et al. . 19981: Ibesch a nd CuzzJ. 
200(j). 



It is clear that in order to constrain the models explaining the crystalline silicates 
in cometary material, accurate knowledge of the abundances of the various com- 
ponents in the dust is required. However, the large differences in the derived crys- 
talline silicate abundances published in the literature so far make progress in this 
area difficult. All studies agree that the abundance of sub-micron sized crystalline 
silicates is very high (although again with considerable spread), but the question 
remains open what the crystallinity of the larger grain population is, where most of 



Email address: mmin@science . uva . nl (M. Min,). 



2 



the mass resides. This is because modeling the optical properties of large, irregu- 
larly shaped particles composed of material with high refractive index is difficult. 



In this study, we combine for the firs t time both the available observations of the 
thermal emission as obtained by ISO ( Crovisier et all ll997l:lLellouch et allll998h 
and th e degree of linear polariz at ion taken fromlGanes 1 et al. ( 1998l):bockers et al. 
( 1999]); lHasegawa etaiT(ll999h : iManset and Bastienl (l200(F and Ijones and G ehrz 
(2000), with a single dust model. In contrast with previous studies we calculate 
the optical properties of the dust grains using a method where we can take into 
account both size and shape effects. With this method, we obtain information on 
the composition and mass of the large grain component. Also we constrain the 
chemical abundances as known from studies of meteorites. Thus, we are able to 
better constrain the dust composition and size distribution in the model than pre- 
vious studies. As we will show, the size distributions of crystalline and amorphous 
silicates are very different, and we are able to rule out the presence of substan- 
tial amounts of large (1-10 pm) crystalline silicates. This results in a much lower 
overall crystallinity than most other studies found. 



In section 2 we explain the fitting procedure and the dust model we employ. The 
results of the best fit model are presented in section 3. The implications of these 
results are discussed in section 4. 



2 Method 



In this section we will outline the method used to interpret the observations of the 
thermal infrared emission and the degree of linear polarization of comet Hale-Bopp 
in terms of the properties of the particles in the coma. The method we use has a few 
important characteristics: 

• Both size and shape effects of the dust grains are taken into account when cal- 
culating the optical properties, i.e. the absorption (emission) cross sections and 
the degree of linear polarization of scattered light. In contrast with previous stud- 
ies, our shape distribution enables us to use the same particle shapes for all dust 
components and grain sizes. This implies that we can, for the first time, include 
the effects of grain size on the optical properties of crystalline silicate grains. 

• The parameters of the size distributions and the abundances of the various com- 
ponents are determined using an 'objective' least squares fitting routine. 

• We require the best fit model to fulfill the constraint that the dust in comet Hale- 
Bopp has the same abundances of the chemical elements of the major solid state 
materials as observed in interplanetary dust particles and meteorites. 
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2. 1 Size and shape of the dust grains 



The optical properties at a particular wavelength, i.e. the absorption and scattering 
cross sections as well as the scattering matrix, of a dust grain are determined by its 
size, shape, orientation, structure and chemical composition. Throughout this paper 
the size of a dust grain is defined as the radius of a volume equivalent sphere, r. We 
use for the size distribution a powerlaw given by 




^min — f — ''max j 



(1) 

elsewhere. 



Here n(r)dr is the number of dust particles of an ensemble with sizes between r and 
r + dr, r mm and r max are the minimum and maximum grain size, respectively, in the 
size distribution; /? is the index of the powerlaw, and C is a normalization constant. 
We choose C so that 



' mas 

I 



4 , 

-7rpr n{r)dr =1, (2) 



where p is the density of the material considered. Hence, n(r)dr is the number of 
dust particles with sizes between r and r + dr per unit mass. The effective radius, 
r eff , and the effective variance, v e ff, are often u sed to characterize a size distribution. 



They are defined by (Mishchenko et al., 2002) 



' max 



r e & = -^r J nr 3 n(r)dr, (3) 



' mm 

' max 



Veff = * 2 f nr 2 (r - r eff ) 2 n(r)dr, (4) 
(G) r 2 eS J 



where 



' max 

= J" nr 2 n 



(G) = nr 2 n{r)dr, (5) 



is the average geometrical shadow of the dust grains. The r eff is simply the surface- 
area weighted mean whereas v e ff is the surface-area weighted variance of the distri- 
bution. 
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The shape of the dust grains is an important parameter determining their opti- 
cal properties. Picture s of interplanetary dust particles (IDPs) show very complex 



shapes and structures (Warre n et all 119941) . To model these complex shapes in de- 



tail requires much computing t ime, which limit s the number of other particle pa 



rameters that can be studied. In lMin et all (|2003) it is shown that when the absorp 



tion cross sections of small particles are considered their shapes can be roughly 
divided into two categories. The first category contains homogeneous spherical par- 
ticles, whereas the second category contains all other particle shapes. Effects of the 
specific particle shape distribution on the absorption properties in the second cate- 
gory are present, but small compared to the differences with homogeneous spheres. 
This implies that the exact shape of the particles is, in a first approximation, not 
important as long as we brake the homogeneity or perfect spherical symmetry of a 
hom ogeneous sphere. This approach is a n example of the statistical approach (see 



e.g. Kahnert et all |2002; Kahnert, 2004). In this approach the absorption proper- 
ties of an ensemble of irregular particles are simulated by the average absorption 
properties of an ensemble of particles with simple shapes. 

In this paper we use a distribution of hollow spherical particles (DHS). In this distri- 
bution a uniform average is taken over the fraction, /, of the total volume occupied 
by the central vacuum inclusion in the range < / < f max , where / max < 1 . All 
particles in this distribution have the same material volume, so that the particles 
with higher values o f f have larg er outer radii. For details regarding this distribu- 
tion we refer to iMin et"all (Eo03 ) . In order to reproduce the wavelength positions 
of the spectral fe atures of crystalline silicate grains, we have to choose f max = 1 



(Min et all 12003). Numerical computations for particles with f = I are not pos- 
sible since these would have an infinitely large outer radius. For particles much 
smaller than the wav elength, integratin g the optical properties up to / = 1 can 



be done analytically dMin et al. . 2003). For larger particles numerical computa- 



tions are necessary. We then choose f max large enough to reach convergence to the 
values for / max = 1. For most cases we considered, it suffices to integrate up to 
/ = 0.98. The optical properties of holl ow spherical particles are calculated us- 



ing a simple extension of Mie theory (Aden and Kerker. 119511) . The calculations 



were done using the l ayer ed sphere code for which the basic ideas are explained in 



Toon and Ac kerman (1981). We would like to stress that we do not allege that hol- 



low spheres are a good approximation for the real shape of cometary dust grains. 
Rather, the properties of the particles in the DHS determining the absorption and 
polarization behavior of the entire ensemble represent in a statistical way those of 
an ensemble of realistically shaped dust grains. 

One of the advantages of the DHS is that the scattering and absorption properties 
can be calculated easily for almost all grain sizes and wavelengths. This means 
that we can calculate the absorption cross sections as well as the degree of linear 
polarization for incident unpolarized light as a function of scattering angle, which 
enables us to fit simultaneously the observed infrared emission spectrum and the 



degree of linear polarization using the same dust parameters. In Mi n et all (|2003) 
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it is shown that good agreement between calculations for the distribution of hollow 
spheres and measurements of the mass absorption coefficients of small forsterite 
grains as a function of wavelength can be obtained. Also, as is shown bv lMin etaL 
(2005]), the degree of linear polarization as calculated using the DHS is in good 
agreement with laboratory measurements of irregularly shaped quartz particles. 



A cometary dust grain consists most likely of a mixture of various components 
such as, for example, olivine, carbon and ironsulfide. In using the model described 
above we assume that the optical properties of such a mixed particle can be rep- 
resented by the average properties of homogeneous particles for each of the sepa- 
rate compone nts. The validity of this approach for core-mantle grains is studied by 
Li et al.l (120021) and shown to be dependent on the shape of the dust grains. In gen- 
eral, whether it is a valid assumption depends on the compactness of the composite 
particle. When the grains are very 'fluffy' aggregates composed of homogeneous 
monomers, the grain components will approximately interact with the incoming 
light as if they were separate. However, in a more compact structure, the effects of 
interaction between the separate components can become visible. 



Throughout the paper we will first average over particle shape (/) which will be 
denoted by (..}. The averaging over particle shape and size will be denoted by <(..)). 



2.2 Thermal emission 



To compute the radiation emitted by an ensemble of monodisperse dust grains in 
random orientation as a function of wavelength we need two ingredients, i) the 
orientation averaged absorption cross section as a function of wavelength, and ii) 
the temperature of the dust grains. When both are known we can calculate the 
observed thermal radiation of a dust grain averaged over all orientations as follows 

nx) = Cabs( ^ (r) . ( 6) 



Here T{X) is the flux density at distance D, A is the wavelength of radiation, T is 
the temperature of the dust grain, B A (T) is the Planck function, and C abs (/0 is the 
orientation averaged absorption cross section of the dust grains at wavelength A. 
The thermal radiation of grains with the temperatures we consider is mainly emitted 
in the infrared part of the spectrum. At these wavelengths the effect of scattering of 
solar radiation on the total observed flux is negligible. 

The temperature of the dust grains is calculated by assuming thermal equilibrium, 
i.e. the energy absorbed is equal to the energy emitted. Note that dust grains with 
equal size r and composition but different shapes and/or orientations can have dif- 
ferent equilibrium temperatures. However, we assume that all dust grains with the 
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same volume and composition have the same temperature determined by the shape 
and orientation averaged absorption properties (as is usually done). 

The coma of Hale-Bopp is optically thin. Therefore, we only need to consider direct 
illumination of the coma grains by the Sun and we may ignore the diffuse radiation 
field caused by, e.g., scattering of solar light by the coma material. We assume the 
Sun to radiate like a black body with a temperature of T = 5777 K. 

Since the absorption cross section, and, therefore, the temperature, depends on the 
size and chemical composition of the particles we determine the temperature for 
each different dust particle size and material separately. Using the size distribution 
given by Eq. (1) the flux per unit mass from dust component j is given by 

r max J 

({fW»; = ^ J nj(r)(C ab!i j(A,r))B,(Tj(r))dr. (7) 

r . minj 

In this equation nj(r) is the size distribution of component j with minimum and 
maximum radii r min j and r maxj ; (c abs j(A, r)^ is the orientation and shape averaged 
absorption cross section of an ensemble of dust grains of material j with size r, 
and Tj(r) is the temperature of the dust grains in this ensemble. Since the coma of 
comet Hale-Bopp is optically thin, the total flux averaged over particle size, shape, 
orientation and composition is simply the sum over the various components 

r model w = J]M j ((nA))) j , (8) 

j 

where Mj is the total mass of dust component j. 



2.3 Degree of linear polarization of scattered light 

In the visible part of the spectrum the radiation from the cometary halo is dominated 
by sunlight scattered once by dust grains. The intensity and polarization of the 
scattered light depend on the angle of scattering and the wavelength. For a comet it 
is possible to obtain measurements of the degree of linear polarization for various 
scattering angles by observing the comet at various moments during its orbit around 
the Sun. All information on the angular dependence of the scattering behavior of 
an optically thin ensemble of dust grains is contained in its 4 x 4 scattering matrix. 
When the size and shape distributions and the abundances of all dust species are 
known, we can calculate the average 4x4 scattering matrix of the ensemble if 
enough data of the refractive index is available. From this matrix we can obtain the 
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degree of linear polarization for incident unpolarized light of dust component j 

((P(a)))j = - } , (9) 

({F n (a)))j 



where a is the phase angle and ({F n k)): is the n, fcth element of the scattering matri x 



averaged over size and shape of the dust grains (for details see Ivan de HulstL ll957). 
It should be noted that for calculations and measurements presented in the literature 
often the scattering angle, 9, is used instead of the phase angle, a. For comets it is 
more convenient to use the phase angle. Since multiple scattering can be neglected 
for comets we have a = 180° - 9. 

In order to calculate the degree of linear polarization of an ensemble of particles not 
only averaged over particle size and shape distributions but also over dust materials, 
we have to average the matrix elements ((F nk ))j. Thus the average polarization is 
given by 

/'model (a) = — — — ■ (10) 

2ujMj{{Fn(a)))j 



2.4 Least squares fitting procedure 



To make a fit to the observations of Hale-Bopp we need to fine-tune the free pa- 
rameters in the model in such a way that we minimize the differences between the 
results of the model computations and the observations in a well defined way. The 
model we constructed has ANj free parameters, where Nj is the number of dust 
species we include in the fitting procedure. By choosing the index of the powerlaw, 
/?, the same for all dust species this reduces to 3N<i + 1 free parameters, namely 



mm,;? ' max,; 5 



fined as 



Mj, j = l..Nd},pj- The error on the infrared emission spectrum is de- 



spec 



N A 



^model(^() _ ^ob&i^-i) 

o- r (Ad 



(11) 



In this equation Ai (i = l..Nj) is a chosen wavelength grid; fobsiAi) is the obser- 
vational value of the flux at wavelength A h and cr T {Ai) is the error of the observed 
flux at wavelength Aj. In order to estimate cr T (Ai) we assume that the error in the 
spectral observations is proportional to the square root of the observed flux. The 
position of the minimum value of x% ec ls independent of the absolute value of the 
error. Note that, since we do not have the values of the absolute errors of the spectral 
measurements, the^ pec is not equal to the reduced ^ and cannot be interpreted 
as the statistical goodness of fit. 
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The error of the degree of linear polarization is defined as 



X 



pol 



EE 



P ohs (ai,Aj) - P mo dei(a:i, Aj) 



(12) 



In this equation P ohs (a h Aj) is the degree of linear polarization observed at phase 
angle a { (i = l..N a ) and wavelength Aj (j = l..N' A ), while a P {ai, Aj) is the error in 
the observed polarization at phase angle or, and wavelength Aj. Note that the number 
of phase angles at which observations are available may vary with wavelength. 

The most straightforward way to define the total x 2 is to take the sum of x% ec and 
X 2 V However, since we have many more measurements of the infrared flux than we 
have observations of the degree of linear polarization and we are unable to compute 
the reduced ;f 2 pec , this would lead to a stronger weight of the spectral measurements 
than the polarization observations. Therefore, we chose to minimize 

(13) 



The minimization of x 2 as defined by Eq. (13) is done using a combination of two 
methods. Since the total spectrum is a linear combination of the separate spectra 
for the different components (see Eq. 8), we are able to separate the fitting problem 
into a non-linear and a linear part. The non-linear part consists of minimizing for 
all r min j, r max j and /?, and the linear part does the minimization for all Mj. 

For the non-li near part of the mini mization we use a genetic optimization algorithm 
called pikaia ( Ch arbonneauL 1 1 9951) . This algorithm tries to find a global maximum 
of an arbitrary function in a large parameter space by using concepts from evolution 
theory. Every set of parameters ({r min j, r max j, j = l..A^},/?) is called an individual, 
the parameters are the 'genes'. The procedure starts with a randomly initialized 
population of N pop individuals and calculates the x 2 for each individual. The indi- 
viduals with the highest values of l/x 2 (lowest values of^- 2 ) are then given the best 
chance to 'reproduce' into the next generation of individuals. This reproduction is 
done by mixing the genes of two individuals into a new individual. When this pro- 
cedure is repeated through several generations, the individuals with low values of 
1 lx 2 will die out and only individuals with high values of 1 /x 2 will survive. In the 
end (after N gen generations) the best individual, which represents the best fit pa- 
rameters, will survive. To ensure convergence also mutations - random variations 
on the para meters - are include d in the algorithm. An extensive description can 



be found in Charbonneau ( 1995). Although the algorithm is not very fast - many 
models need to be calculated - it is very robust in the sense that it will (almost) 
always find the global maximum, whereas other optimization codes frequently end 
up with a local maximum. 
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Genetic algorithm 



First generation 
(randomly initialized) 




Given the set of parameters 

({'■„,,„ j. <-„,„j..; = i-JV} 




[Vj, j = 1,^1 are calculated by minimizing 
Xtpa: us i n S a linear least square fitting 
procedure 



Using the (|r m ,„ j, rmmJ , j = l..N],fi) 
and the obtained values for {Vj,j = l,N] 
XL is com P uted 




Fig. 1 . A schematic representation of the fitting procedure. 

For every individual in the genetic algorithm we determine the best values for the 
Mj by minimizing the x% e c using a linear least squares fitting procedure. Linear 
least squares fitting amounts to solving an overdetermined matrix equation in a 
least squares sense. To ensure that all Mj are positive we need a robust linear least 
squares fitting procedure with extra linear equality and inequality constraints. We 
use the subroutine dlsei from the slatec library 1 . The inequality constraints are 
used to ensure that Mj > 0. The equality constraints are used to constrain the chem- 
ical abundances as will be explained in section 2.6. Using the Mj thus derived, we 
calculate ^p Gl and^ 2 . Although in this way we might miss the absolute minimum, 
we can use this approach since the infrared spectrum is more sensitive to the exact 
dust composition than the degree of linear polarization is, which is more sensitive 
to the size distribution. 

In order to estimate the errors on the derived abundances of the dust species we 
consider all individuals from all generations with a value of x 2 smaller than 1 . 1 
times the minimum value of^- 2 . The value 1.1 was chosen such that all fits within 
this range are still in reasonable agreement with the observations. 

The entire fit procedure is schematically outlined in Fig. 1. 



The slatec library is publicly available for download at: http://www.netlib.org/slatec/ 
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2.5 Dust components 



In this section we will discuss the dust components we chose for our model com- 
putations. Based on previous studies and on the composition of IDPs we make the 
following selection of components. 

• Amorphous Carbon (C). We assume that most of the Carbon atoms that are in 
the solid state phase will be present as amorphous carbon. The emission spec- 
trum of amorphous carbon gives a smooth continuum contribution without clear 
spectral structure. From the in situ measurements by the Vega spacecraft when it 
encountered the coma of comet Halley, we know th at approxima tely half of the 
total available C is present in the solid state phase (iGeissL 119871) . Carbon prob- 
ably acts as a matrix in which the other materi als are embedded. The refractive 
index as a function of wavelength is taken from lPreibisch et al.l (11993). 

• Amorphous Olivine (Mg2 A -Fe 2 -2A-SiC)4). Amorphous olivine is one o f the most 



abun dant dust species in circumstellar and interstellar material (see e.g . Molster et al 
1999; Bo uwman et alll2001uLi and DraineL 12001: For rest et al.L l2004: Kempe r et al. 



2004). The sili cate dust in the diffuse interstellar medium is dominated by amor- 



phous olivine (IKemper et all 12004). The emission by small amorphous olivine 
grains shows broad spectral features at 10 and 20 jjm. Both are detecte d in the 



spectru m of Hale-Bopp. T he data for the refractive indices are taken from Do rschner et al 
dl995h . 

Amorphous Pyroxene (Mg^Fei^SiOs). Amorphous pyroxene is spectroscopi- 
cally hard to distinguish from amorphous olivine. Pyroxene is often found in 
IDPs. The emission by small amorphous pyroxene grains shows a spectral struc- 
ture which is similar to that of amorphous olivine grains. However, the maximum 
of the 10 /im feature is shifted towards slightly shorter wavelengths and the shape 
of the 20 /im feature is slightly different. The r efractive indices for amorphous 
pyroxene are taken from Dorschner et al.l (11995b . 

Crystalline Forsterite (Mg 2 Si0 4 ). Crystalline forsterite is the magnesium rich 
end member of the crystalline olivine family. Experiments indicate that when 
amorphous olivine is annealed under certain conditions, the iron is removed 
from the lattice structure and crystalline forsterite is formed. The most impor- 
tant resonances in the emission spectrum of small crystalline forsterite grains are 
located at wavelengths 1 1 .3, 19.5, 23.6 and 33.6 /urn. These features are observed 
in the spectra of, for example, AGB star s and protoplanetary disks (see e.g. 
Waters et al. , 1996; Bouwman et al., 2001) and are clearly visible in the spec- 
trum of Hale-Bop p dCrovisier et allll997l). F or the refractive indices of forsterite 
we use the data of IServoin and Pirioul ( 1973 ). 

Crystalline Enstatite (MgSi0 3 ). Crystalline enstatite is the magnesium rich end 
member of the crystalline pyroxene family. Under the right conditions, enstatite 
can form from a reaction between forsterite and silica. Also, from IDPs there 
are indicat ions that some enst atite grains formed directly from gas phase con- 



densation (Bradle y et all 119831) . The spectrum of small enstatite grains shows 
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prominent features at 9.3, 10.5, 19.5, 28 and 44.5 //m. Although we have no di- 
rect evidence for the presence of most of these features, we include enstatite in 
the fitting procedure because enstatite features have b een reported in the g round 



based spectrum of Hale-Bopp when it was at 1.2 AU (Wooden et al., 1999). The 
data for the refractive indices are taken from ljager et aL ( 19981) . 



• Amorphous Silica (Si02). Laboratory measurements show that when amorphous 
silicates are an nealed to form crys talline silicates, a certain amount of silica is 
also produced (IFabian et all 12000). The emission spectrum of small amorphous 
silica grains shows features at 9, 12.5 and 21 pim. We take the refracti ve indices 
as a function of wavelength measured by Spitzer and Kleinma^ dl960l) . 

• Metallic Iron (Fe). When amorphous silicates are annealed to form magnesium 
rich crystalline silicates, the iron is removed from the lattice. The spectral sig- 
nature of metallic iron grains is very smooth but slightly different than that of 
amorphous carbon. Although it is not crucial in obtaining a reasonable fit, we 
include the possibility of metallic iron grains in the fitting proce dure for com 



pleten ess. The refractive indices of metallic iron are taken from Hennin g et al. 
dl996h . 



• Iron Sulfide (FeS). In IDPs all available sulfur is present as iron sulfide. The 
spectral structure of emission from iron sulfide grains shows a broad feature 
around 23 jjm. The fact that this feature is not prominent in the spectrum of Hale- 
Bopp is probably due to the fact that the iron sulfide grains are relatively large, 
which reduce s the spectral structure significantly. We use the refractive indices 



measured by Begem ann et al.l (11994) 



To decrease the number of free parameters, we have chosen equal size distributions 
for both amorphous silicate species (amorphous olivine and amorphous pyroxene) 
and equal size distributions for both crystalline silicate species (forsterite and en- 
statite). Furthermore ft, the index of the powerlaw, is the same for all materials. 



2.6 Chemical abundance constraints 



The number of free parameters in the method as described above is very large. 
When we apply no extra constraints we encounter a large set of solutions that all 
have a more or less equal but with totally different size distributions and ma- 
terial abundances. We also run the risk of obtaining a best fit solution with very 
implausible values for the fit parameters. To avoid these problems we constrain the 
solution by requiring that the major elements for solid state materials (Si, Mg, Fe 
and S) are all in the solid phase and that their chemical abundances are the same 
as those found in meteorites. Furthermore, we take the abundance of Carbon in the 
solid phase to be half the solar abundanc e, cons i stent w ith in situ measurements 
of the dust in the coma of comet Halley (Geiss. 119871) . This gives us four extra 
constraints on the model. These extra constraints prove to be sufficient to obtain 
consistent results with the fitting procedure. 
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Table 1 



V llV^llllV^lLl V^lV^lllV^llL 




Solar 


f^nn sir flints 


C/ol 




y.jj 


4.0 / 


Mg/Si 


1.05 


1.07 


1.05 


Fe/Si 


0.87 


0.89 


0.87 


S/Si 


0.44 


0.60 


0.44 



Abundance constraints as a pplied in the fitting procedure. These values are taken from 
Grevesse and Sauval dl998h . For Carbon we take half of the s olar abundan ce, consistent 
with in situ measurements of dust in the coma of comet Halley ( GeissL 1987 ). 



We constrain the abundances of C, Mg, Fe and S relative to Si. For the Carbon 
abundance we take half of the total solar abundance to be in the solid phase, 
the remaining Carbon is in the the gas phase. Magnesium, iron, sulfur and sili- 
con are ass umed to be complete l y in th e solid phase. The abundances we use are 
taken from lGrevesse and Sauval dl998l) and are summarized in Table 1. Note that 
the sulfur abundance measured in meteorites is lower than the solar abundance. 
Thus by taking the constraints from meteorites for the solid state particles we as- 
sume there is also sulfur in the gas phase . This is consistent with observations 
(Bockele e-Morvan et al. . 2000l: Irvine et al. . 2000l) . All abundance constraints are 
incorporated as linear equality constraints to the linear least squares fitting part of 
the minimization procedure. 



To use the chemical abundance constraints we have to introduce an extra free pa- 
rameter, x, the magnesium fraction in the amorphous silicates. In section 3.2 we 
will show that the bulk of the material consists of olivine, pyroxene, carbon and iron 
sulfide. Using this information we can already make a simple but reliable estimate 
of the value of x by adopting these four species only and applying the abundance 
constraints discussed above (as this implies that we have four constraints and four 
unknown parameters the material abundances are uniquely defined). The results of 
this simple calculation are summarized in Table 2, and give x = 0.7. Note that 
here it is not possible to distinguish between amorphous and crystalline material. 
The results from the fitting procedure will be slightly different due to a different 
composition of crystalline and amorphous silicates and the fact that in the fitting 
procedure we also added silica and metallic iron. 
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Dust component 


Chemical formula 


Abundance (Mass %) 


Carbon 




z4.o 


Olivine 


Mg2xFe 2 -2xSi0 4 


33.4 


Pyroxene 


Mg Y Fe 1 _ t Si0 3 


25.0 


Iron Sulfide 


FeS 


17.0 



Table 2 

Abundances as calculated from chemical abundances found in meteorites. In order to satisfy 
the chemical abundances given in Table 1 , we have to take x = 0.7. 

3 Results 



3.1 Observations 



The spectroscopic observations we use are the infrared spectra obtained by the 
Short Wavelength Spectrometer (SWS) and the Lon g Wavelength Spectrometer 



(LWS ) on board the Infrared Space Observatory (ISO) (Cr ovisier e t all 1997t lLellouch et al. 



19981) . These spectra were taken when the comet was at 2.9 AU distance from the 



Sun and 3 AU distance from the Earth. The SWS and LWS spectra have a small 
overlapping wavelength range (42 /im < A < 45 /im). Since the LWS has a larger 
beam size, it catches emission from a larger part of the coma, resulting in a higher 
absolute flux level. We assume that the properties of the dust causing the emission 
does not change as a function of distance from the core of the comet, so we can 
simply scale the LWS spectrum to match the SWS spectrum in the overlapping 
wavelength region. We have to note here that there are indications that the size dis- 
tribution or the compa ctness of the particles v aries slightly as a function of position 



in the coma (see e.g. Kolok olova et all 120041 and references therein). When the 
particles move away from the comet nucleus, the particles might fall apart resulting 
in more fluffy, or smaller structures. Therefore, one might argue that taking into 
account the extended region covered by the LWS might bias our results towards a 
slightly higher fraction of small grains. However, we believe that this effect is only 
minor since these differences are largest when considering the region very close 
to the coma, which is covered by both the SWS and the LWS (Kolokol ova et all 



2004). For the fitting procedure we used the wavelength range from 7 to 120 fim. 



For the observational data of the degree of linear polarization we used the com- 
bined measurements from various studies in the optical to near infrared part of the 



spectr u m. The measurements were taken from Ganesh et al. ( 1998); Jockers et al 
( 1999); Hasega waetaD dl999h : iManset a nd Bas tienl d2000h and Ijones an d Gehrz 



(2000). The observations also provide the errors o>. Combining the observations 
we have polarization data at twelve different wavelengths for various phase angles. 
We note that in order to obtain observations at different phase angles, the comet 
has to be observed at different phases during its orbit around the Sun. Therefore, in 
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Fig. 2. Infrared spectral energy distributions of the best fit model (black line) together with 
the SWS and the LWS observations (gray line). 

order to model all these observations using a single dust model, we have to assume 
that the composition and size distribution of the dust is more or less constant at 
different phases. Since we consider the degree of linear polarization of scattered 
light, variations in the total dust mass in the coma are not important. 

3.2 Best fit model 

We minimized^ 2 as defined by Eq. (13) using the method described in section 2.4. 
The Ai were chosen on a logarithmic grid. Throughout the fitting procedure we fix 
the value of f max for each dust species. From the positions of the crystalline silicate 
resonances we already know that for these materials we have to choose the most 
extreme shape distribution parameters (/ max = 1). The other materials are chosen to 
have equal values of / max . For these materials, we obtain an optimum value / max = 
0.8. The parameters are shown in Table 3. The infrared spectrum corresponding to 
the best fit model is shown in Fig. 2 together with the measurements. The resulting 
curves for the degree of linear polarization as functions of the phase angle, together 
with the observations, are shown in Fig. 3. The emission spectra of the separate 
dust components are plotted in Fig. 4. 

Figs. 2 and 3 show that the observations of both the infrared spectrum and the de- 
gree of linear polarization can be reproduced remarkably well using the same dust 
model. In other words, we do not need quite different models for the infrared and 
optical parts of the spectrum. The differences between the observed and predicted 
infrared spectrum are most probably mainly caused by uncertainties in the refrac- 
tive index data, and the assumptions on which the model is based. For the degree 
of linear polarization we notice that the negative polarization branch at A = 365 nm 
(and to a smaller extend at 484.5 and 684 nm) at small phase angles is not repro- 
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Material 


Chemical Formula 


Bulk density 


fmax 


Volume 


(min / max) 


Mass 


^rnin 


r max 


P 






(g/cm 3 ) 






(%) 


(%) 




Gum) 




Amorphous Carbon 


C 


1.80 


0.8 


39.1 


(39.0/39.3) 


23.8 


7.0 


11.6 


-3.48 


Amorphous Olivine 


Mg 2 ^Fe 2 -2^Si04 


3.71 


0.8 


25.7 


(22.9/32.7) 


32.3 


0.01 


92.8 


-3.48 


Amorphous Pyroxene 


Mg,Fei_,Si0 3 


3.20 


0.8 


18.3 


(7.8 /23.5) 


19.8 


0.01 


92.8 


-3.48 


Amorphous Silica 


Si0 2 


2.60 


0.8 


2.2 


(0.8 / 5.6) 


1.9 


6.5 


7.5 


-3.48 


Crystalline Forsterite 


Mg 2 Si0 4 


3.33 


1.0 


2.7 


(2.5 / 3.0) 


3.1 


0.05 


0.1 


-3.48 


Crystalline Enstatite 


MgSi0 4 


2.80 


1.0 


1.3 


(0.9 / 1.7 ) 


1.2 


0.05 


0.1 


-3.48 


Metallic Iron 


Fe 


7.87 


0.8 


0.6 


(0.3 / 1.1) 


1.5 


0.2 


0.3 


-3.48 


Iron Sulfide 


FeS 


4.83 


0.8 


10.1 


(10.0/10.1) 


16.4 


0.8 


4.3 


-3.48 



Table 3. Results for the best fit model. The error estimates are obtained by considering all possible fits with^ 2 < 1 . ly ^, in ■ The total dust mass in 
beam of the SWS as found from the fit is 4.6 • 10 9 kg. 
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Fig. 3. The degree of linear polarization from the best fit model as a function of phase angle, 
a, in degrees (solid curves) based on both the spectral and the polarization measurements. 
The measurements o f the linear po l arizati o n are indicated by dots a nd were taken from 



Ganesh et al. ( 1998); Jockers et al 


] <ES99); 


Haseaawa et al. ( 1999); Manset and Bastien 


fcOOC 


) and Jones and Gehrz 


(2000 





duced satisfactorily by the model. This is most likely due to the spherical symmetry 
of the shapes we employed. In addition a discrepancy between the model and the 
observations occurs at A = 1250 and 1650 nm, which is possibly also connected 
to the adopted spherical particle shapes. The scattered light is dominated by the 
contributions from the silicate (olivine and pyroxene) and ironsulfide grains. The 
contribution to the scattered light from carbon grains is negligible. Although the 
abundance of silicates is much higher than that of ironsulfide, the scattering caused 
by ironsulfide grains is comparable to that of the silicates due to its high scatter- 
ing efficiency. Especially at near infrared wavelengths (the A = 1650 and 2200 nm 
measurements) the scattering is dominated by ironsulfide grains. 

The best fit model we present combines for the first time observations of the SWS 
infrared spectrum (7 - 44jum), the LWS infrared spectrum (44 - 120/^m) and the 
degree of linear polarization at several wavelengths in the optical to near infrared. 
This results in a better constrained dust model. The spectral structure of the ther- 
mal emission in the SWS part of the spectrum provides crucial information on the 
composition of the dust. The LWS part of the spectrum combined with the degree 
of linear polarization provides information on the size of the dust grains. We also 
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Fig. 4. The contributions of all small (dots), and large (dashes) grain components to the total 
spectrum (upper panel). The gray line indicates the observations, the black line indicates 
the best fit model. Also plotted are the contributions to the emission spectrum of the various 
dust components (lower panels). 

constructed fits excluding some of the observations from the model. An attempt to 
fit only the SWS part of the spectrum resulted in an underestimate of the fraction of 
large grains, which in turn resulted in a hi gher fraction of crystalli ne silicates, more 
comparable to that found in, for example. Bouwrnan et alJ (|2003). Fitting both the 
SWS and the LWS part of the spectrum without the linear polarization resulted in 
a dust composition only slightly different (i.e. within the given error bars) from 
that presented in Table 3. The size distribution is affected more significantly. Also, 
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using observations of the linear polarization at less wavelength points changed the 
parameters of the best fit model. For example, the best fit model when the observa- 
tion of the degree of linear polarization at A = 1650 and 2200 nm are removed from 
the model contains no amorphous pyroxene and 8.5% amorphous silica. This can 
be e xplained in part by the interchangeability of olivine and pyroxene grains (see 
also lvan Boekel et all 120051) . The size distribution of this best fit model is changed 
such that the upper size limit of the amorphous olivine component is only 30 /im 
but the slope of the size distribution is more weigthed towards the larger grains, 
a = 3.0. The other parameters are only affected mildly. 



The total dust mass in the SWS beam derived from t he best fit model is 4.6 • 10 9 kg. 
This dust mass is comparable to that estimated by Bouwrnan et all (120031) which 
is 4.2 • 10 9 kg. Note that the dust mass derived in this way is a lower limit on 
the real mass of the solid state material in the coma of Hale-Bopp since the mass 
most likely resides predominantly in the very large grains. Using data obtained at 
submil limeter wavelengths, whic h provides information on millimeter sized dust 
grains, Jewitt and Matthews! d 19991) derive a total dust mass of ~2 • 10 11 kg within a 
beam size comparable to that of the ISO SWS. 



In our model the temperature of all dust species is determined self-consistently. 
However, from the ratio of the strengths of the forsterite features it can be seen 
that the temperature of the forsterite grains is higher than would b e determined 
from thermal equilibrium calculations using pure forsterite grains ( Harker et all 
iooi . Pure forsterite grains are not very efficient absorbers in the UV and the op- 



tical part of the spectrum where a large part of the solar energy is emitted. They 
are very efficient emitters in the infrared. Therefore, pure forsterite grains will be 
cold compared to other dust species. The fact that they are observed to be rel- 
atively warm is probably an e ffect of thermal contact between the various dust 
species (IBouwman et alll2003r) . An aggregated structure where all dust species are 
in thermal contact is also consistent with pictures of interplanetary dust particles. 
To calculate the optical properties of an aggregated structure in a completely con- 
sistent way is very computationally intensive which makes it extremely difficult 
to examine the entire parameter space of dust abunda nces and grains sizes. Some 
work on this has been done by [Moreno et all d2003l) considering fixed values of 
the material abundances. To simulate thermal contact we polluted the forsterite and 
enstatite grains with 3% of small metallic iron inclusions. We calculated an effec- 
tive refractive index using t he Maxwell- Garnet effective medium theory (see e.g. 



Bohren and Huffman, 1983). When we pollute the crystalline silicates in this way, 



the equilibrium temperature is in agreement with the observations. 
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3.3 Dust composition and size distribution 



When discussing the size distribution of the dust we have to consider the grain sizes 
that our analysis is sensitive to. For grains smaller than a few micron (< 3 //m) we 
have a strong spectroscopic diagnostic. At relatively short wavelengths (around 
~10/jm) we are mainly sensitive to the composition of these small grains. If we 
go to longer wavelengths, we are also sensitive to the composition of larger grains. 
Using the SWS and LWS range from A = 7 - 120/im, we have a spectroscopic 
diagnostic for the composition of grains with a volume equiva lent radius up t o 



~10-15//m. Although larger grains do show spectral structure (Min et all 12004). 
their emission efficiency is too low to be detected. 

When we compare our best fit parameters with those obtained by others ferucato et all 



1999; 



2002; 



Gald emard etat[ll999l:IWooden et all[l999l:lHavward et alll200(ilHarker et al 
Bouwman et all 2003 ) there are a few differences. First of all the amount of 



crystalline silicates is much smaller than that found in most of these studies. This 
is caused by the fact that previous studies mainly considered the small (submicron 
sized) grain component. In our model, we find that in order to reproduce the spec- 
tral features, the crystalline silicate grains have to be very small and thus they have 
a high emission efficiency. From, i.e. the LWS spectrum, we find that the amor- 
phous grains are relatively large, and thus emit less efficiently. Thus, if one only 
considers the small grain component, the crystallinity is increased with respect to 
our findings. The crystalline grains have to be small in order to reproduce the ratios 
of the strengths of the different emission features. However, these ratios are also 
influenced by the temperature of the dust grains. Fortunately, we can distinguish 
temperature effects from grain size effects in the following way. In determining 
the absorption cross sections of dust particles the most important parameter is \mx\ 
where m is the complex refractive index and x = 2nr/A is the size parameter of the 
dust grain. When \mx\ « 1 the grains are in the Rayleigh domain and strong spec- 
tral emission resonances occur. When \mx\ increases, the spectral features decrease 
in strength. For very l arge grains (\mx\ » 1) the emission features will change 
into emission dips (see lMin et al. . 2004h . When the grain size is increased, first the 
features caused by resonances with high values of \m/A\ become weaker and then 
the features with smaller values of \m/A\. This means that this effect depends on the 
wavelength but also on the refractive index; the strongest features (with the highest 
values of |m|) will go down first. However, when going from high temperatures to 
low temperatures, the effect on the strength of the features shifts with the maximum 
of the underlying Planck function, so this effect only depends on the wavelength 
position of the features. This difference allows us to distinguish between size and 
temperature effects by carefully looking at the feature strength ratios. 



To test if the ratios of the forsterite resonances can be explained using temperature 
effects instead of a difference in grain size between amorphous and crystalline sili- 
cates, we calculated the strengths of the 19.5,23.6 and 33.6/mi features above the 
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•^19.5/ ^23.6 

Fig. 5. The ratios of the peak strengths at 19.5, 23.6 and 33.6 fxm of forsterite calculated for 
various grain sizes and temperatures. The star indicates the measured peak strength ratios 
of comet Hale-Bopp. Grains smaller than 0.1 jum are in the Rayleigh limit. In this limit, the 
shape of the spectrum is independent of the particle size. 

local continuum using the DHS for different grain sizes and temperatures. These 
strengths are denoted by S 19.5, S 23.6 an d 5*33.6 respectively. The strength ratios are 
plotted in Fig. 5 for various grain sizes and temperatures. The temperatures were 
varied from 10 - 1500K. Also plotted in Fig. 5 is the peak ratio measured in the 
spectrum of Hale-Bopp. We see that there is no corresponding set of temperature 
and grain size in Fig. 5 that reproduces the ratios measured in Hale-Bopp. To get 
the best fit we need to go to very small and relatively warm forsterite grains. The 
fact that the forsterite has to be warm also indicates that the grains must be small 
and most probably are in thermal contact with (at least) a strong absorbing material 
like carbon or iron. 

The fact that we cannot find a set of temperature and grain size with peak ratios in 
Fig. 5 corresponding to those observed in the infrared spectrum could be caused by 
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Fig. 6. Infrared spectral energy distributions of the best possible fit when fixing the maxi- 
mum volume equivalent radius, r max cryst , of the crystalline component to 5 /mi (black line). 
This fit to the SWS and LWS observations (gray line) therefore does not represent our best 
fit model (given in Fig. 2, which uses r max cryst ~ 0.1 /urn), but is of a much poorer quality. It 
is clear that large crystalline silicate grains cannot be present as they cannot reproduce the 
observed peak ratios of the resonances. 

the effect that in the spectrum of Hale-Bopp the 33. 6 //m feature is partly blended 
with an enstatite feature resulting in a weaker feature. If we would take this effect 
into account, the point would shift downward in Fig. 5. Another explanation could 
be that the forsterite in Hale-Bopp is slightly contamin ated with iron, which results 
in a slightly weaker 33.6 jxm feature (Ko ike et allll993h . 



As a test we made a fit to the SWS and LWS spectra fixing the upper size of 
the crystalline silicate grains. We have tried to make a fit to the spectrum using 
Vax,cryst = 2,5 and 10/mi, respectively, and varying all other parameters. The re- 
sulting model spectra did not satisfactorily reproduce the measured ISO spectra. In 
order to obtain a reasonable fit to the observed forsterite features, we had to employ 
a power law for the size distribution more biased towards small grains, fi = -3.6. 
However, in all these model fits the strength ratios of the forsterite features were 
poorly reproduced. In Fig. 6 we plot the resulting best fit model for the case when 
Vax,cryst = 5 yum. We note that the emission spectrum is not sensitive to the very 
large forsterite grains. Although these grains still display significant spectral struc- 



ture (Mi n et all l2004f) . their emission efficiency is small. Therefore, it is possible 



to have a bimodal size distribution of crystalline silicates, in which only very small 
and very large crystals are present while the intermediate sized grains are absent, 
and still reproduce the infrared spectrum of Hale-Bopp. However, such a size dis- 
tribution is very unlikely. 



When comparing our results to those obtained by other studies, we have to be care- 
ful with the definition of crystallinity. Here, we mean by crystallinity the mass 
fraction of crystalline silicates compared to the total dust mass. In previous studies, 
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sometimes the crystallinity is defined as the mass fraction of crystalline silicates 
compared to the total mass in the silicate component. In our fit, approximately 40% 
of the dust mass is contained in non silicate materials like carbon and ironsulfide. 
When comparing with other studies we have to correct for this. 

Since previous studies of the mineralogy of the dust in the coma of Hale-Bopp 
focused on the SWS part of the spectrum, little information was available on the 
large grain component. As can be seen from Fig. 4 the LWS spectrum (long ward of 
~44 /urn) displays a long wavelength slope that can only be explained using large 
grains at approximately blackbody temperature (see also lLellouch et al. i ll998h . We 



have just shown that this large grain component cannot be crystalline. Since the 
large grains contain most of the mass, the resulting abundance of crystalline sil- 
icates is small, 4.3% (Table 3). We should note here that part of the crystalline 
silicates are probably not spectroscopically detectable since these are hidden inside 
large amorphous silicate or carbon grains. Therefore the 4.3% crystallinity we de- 
rive is a lower limit. We can make a quantitative estimate of the true fraction of 
crystalline grains by only considering grains smaller than a certain size, for which 
we assume that we observe the properties of the entire grain. In order to assure 
that the resulting mixture still has solar abundances, we fix the total abundance of 
silicates (crystalline and amorphous olivine and pyroxene). Considering only the 
silicate component, we compute the crystalline over amorphous ratio excluding the 
larger silicate grains. We then use this crystalline over amorphous ratio to compute 
the fraction of crystalline silicates in the total mixture. We choose a volume equiv- 
alent radius of 10 fim as the maximum grain size, which represents approximately 
half of the total dust mass. The crystallinity derived in this way depends on the as- 
sumption of this grain size, and detailed calculations are required to obtain a better 
estimate of the appropriate size. We thus derive a crystallinity of ~7.5% (~5% of 
crystalline forsterite and ~2.5% of crystalline enstatite). Note that most previous 
studies on the SWS spectrum of Hale-Bopp only included forsterite, no enstatite. 
In contrast with our findings, these studies typically find forsterite fractions of ~20- 
30%. The fraction of crystalline silicates in our best fit model increases when we 
consider only smaller grains. The submicron grain component contains approxi- 
mately 30% forsterite and 14% enstatite. This fraction of forsterite is consistent 
with the resu lts from previous stu dies that focused on the submicron grain compo- 
nent (see e.g. Brucato et al. . 1999b . 



It can be noted from Table 3 that the dust components with a high abundance reside 
in relatively large grains. Although the large grains contribute much to the total 
dust mass, they contribute only little to the infrared emission at short wavelengths. 
Therefore, we conclude that the long wavelength part of the spectrum and the de- 
gree of linear polarization contain crucial information on the large dust grains, and, 
therefore, on a large fraction of the total dust mass. 

We can calculate the effective radius and variance, r e g and v e ff, of the dust grains 
in our model using Eqs. (3) and (4). For our best fit model r eff = 1.0 //m and 
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v e ff = 16.8. Typical values for the effective volume equ ivalent radii computed from 
the size distributions found i n other studies are 0.5 /mi dBouwman et al. . 2003b and 



0.6 /mi (IMoreno et all 120 03). The larger r eff that we find is most likely due to the 
additional information contai ned in the long waveleng th slope. Employing a simple 



model for porous dust grains JLi and Gree nberg ( 1998) argue that the average grain 



size in Hale-Bopp might be much larger (~8-25 //m) and that the fraction of sub- 
micron grains might be negligible. Perhaps cometary dust grains are indeed very 
fluffy aggregates of small particles. However, it is probably difficult to discrimi- 
nate between a large fluffy grain and a large number of small separate grains using 
spectroscopy or observations of the degree of linear polarization. 

The fact that the most abundant dust species reside in the largest grains is a natural 
consequence of the formation mechanism of the larger dust grains. Larger dust 
grains are believed to form by coagulation of smaller grains. The materials that are 
more abundant, can form large clusters of this material, whereas the less abundant 
species may be distributed in the aggregate as separate monomers. If the aggregate 
is very fluffy, the low abundance of the latter species causes the distance between 
the inclusions to be on average relatively large. Thus, one could argue that these 
inclusions will not interact very strongly with each other, and their resulting optical 
properties will be those of small particles. On the other hand, the larger clusters 
formed by the more abundant dust species will produce optical properties more like 
those of large grains. This hypothesis may be tested by simulations of the optical 
properties of fluffy composite particles with extreme abundance differences. 

We conclude that the crystalline silicates in the coma of comet Hale-Bopp are sub- 
micron sized. This is much smaller than the typical grain size found for the other 
dust components. In agreement with these fi ndings, the crystallin e silicates found 
in IDPs are predominantly submicron sized (Bra dley et all 1 1999). This difference 
in grain size between the amorphous and crystalline components results in a rela- 
tively low overall abundance of crystalline silicates. The abundance of crystalline 
silicates and the difference in grain size are important constraints when considering 
models of the processing and dynamical history of dust in the early solar system. 



4 Discussion: Origin and evolution of cometary dust 



Since comets are small bodies that have been frozen all of their lifetime, the comet 
material is expected to have undergone little processing since the comet's forma- 
tion. In larger bodies, like planets, the material has experienced severe parent body 
processing erasing the information on the original dust grains. Therefore, comets 
contain a unique diagnostic of the mineralogy of the dust in the protoplanetary disk 
at times when the planets and comets were formed. The fract ion of crystall i ne sil- 
icates in the diffuse interstellar medium (ISM) is very low. iLi and Draind (|2001) 
derive an upper limit for the crystallinity of the dust in the ISM of 5%, while 
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Kem per et al. (2004) even derive a n upper limit of 0.4%. In addition, the crys 



tallinity of the ISM dust derived by van Boekel et al. (2005) is ~1%. This implies 



that the dust in Hale-Bopp has undergone at least some processing in the solar neb- 
ula, before it got incorporated into the comet. In order to form crystalline silicates, 
the amorphous silicates have to be heated above the glass temperature of ~1000 K. 
Comets form in the outer regions of the disk where the temperatures are low enough 
for water ice to exist < 160K. There are two possible ways to explain the pres- 
ence of crystalline silicates in cometary dust grains. The first is that the amorphous 
silicates are crystallized by thermal annealing in the hot inner regions of the pro- 
toplanetary disk. The crystalline silicates are then mixed out to the regions where 
the comets form. Another possibility is that the crystalline silicates are produced 
local ly. Processes that have b een proposed fo r this are, for example, shock anneal- 
ing dHarker and DeschL 120021) and lightning (Pilipp et all 1 19981: IDesch and Cuzzl 
2000 ). Recent evidence for radial mixing in protoplanetary disks is presented in 
Ivan Boekel et all (120041) on the basis of interferometric measurements of the dust 
in the inner disk regions. In that paper it is found that the crystallinity in the in- 
ner disk regions is higher than that in the outer disk regions. This difference, along 
with a varying forste rite over en statite ratio, is consistent with predictions from 
radial mixing models dGailL 12004). 

Several studies have tried to explain the presence of crystalline silicates in comets 
by considering radial mixing. In these studies only the evolution of the silicate com- 
ponent is computed. For comparison with our computations, we therefore have to 
correct for the additional components included in the fit. The fraction of crystalline 
silicates in the silicate compone nt (amorp hous and crystalline olivine and pyrox- 



l ent (i 
l Eail 



ene) is approximately 1 2.5% . In Gail ( 2004) a detailed model for the protosolar disk 



is presented in which dust chemistry, thermal annealing and radial mixing are in- 
corporated. The computations in this paper are for a stationary model, so an equilib- 
rium situation is calculated. The crystallinity that follows from this detailed model 
depends strongly on the distance to the star. Close to the star, the crystallinity is very 
high (~70% at 3 AU), while in the comet forming region (~20 AU) the crystalline 
fraction equals approximately 24% (8% crystalline olivine and 16% crystalline py- 
roxene). The author concludes that the results from the stationary model cannot be 
extrapolated to distances beyond 20 AU, and can only be considered an approxi- 
mation for the mineral co mposition of the inner 20 AU at several times 10 5 years. 



Wehrstedt and Gail (2002) present a time depend model of radial mixing. From the 
results of t his paper, it is apparent that the equilibrium situation as computed in 
Gail (2004) is already attained after ~10 5 years. This is slightly shorter than the an- 
ticipated typical time scale for comet formation (which i s a few times 10 5 ye a rs, see 



IWei denschillingl fl997 ). The crystallinity computed by We hrstedt and Gai3 d200 



at a d istance of 30 AU is approximately 7% after 10 6 years . iBockelee-Morvan et al 
(2002) also present a model describing time-dependent radial mixing in a proto 



planetary disk. In this paper, three different solar nebula models are presented, a 
warm, a nominal, and a cold model referring to the temperature structure in the so- 
lar nebula. This temperature structure is set by the viscosity parameter a; the higher 
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the values of a the lower the temperature. They derive an extremely well mixed 
nebula, in which, after ~10 6 years, the crystallinity at distances > 10 AU is inde- 
pendent of the distance. They arrive at a final crystallinity in the outer solar system 
(> 10AU) of approximately 58, 12 and 2% according to the warm, nominal, and 



cold s olar nebula mod el. Note that the parameters chosen by Wehrstedt and Gail 



(2002) and Gail d2004 corre spond to the warm to nominal solar nebula model of 



Bockelee-M orvan et al.l d2002l) . The predict ed crystallinity in th e comet forming 



region from both Gail (2004) and B ockelee-Morvan et al.l (|2002) shows that ther 



mal annealing and radial mixing are more than sufficiently efficient mechanisms to 
explain the crystalline silicates in Hale-Bopp. The most important free parameter 
in the above models is the viscosity parameter. The crystallinity of cometary dust 
can be used to constrain this parameter, providing a better insight in the dynamics 
of protoplanetary disks. While the crystallinity derived for Hale-Bopp in previous 
studies could only be explained employing a viscosity parameter that is represen- 
tative for a warm solar nebula model, the crystallinity we derive for Hale-Bopp is 
consistent with a viscosity parameter that is typical for an approximately nominal 
solar nebula and a formation distance of some 30 AU from the central star. 

The crystalline silicates we find are all submicron sized, consistent with studies 
of IDPs, in which th e crystalline inclusions are predominantly submicron sized 
(Bradley etafl[l999). There are two possible explanations for this. The first ex- 



planation is that the crystalline silicates are formed before efficient grain growth 
sets in. Due to the low crystallinity this would result in a dust grain that has the 
crystalline silicates scattered in the aggregate as small separate inclusions. How- 
ev er, this possibility can likely be excluded on the basis of the results presented 



by |van Boek el et al.l (120051) . From the analysis of a large sample of protoplanetary 
disks, it is concluded that grain growth occurs before efficient crystallization sets 
in. Another explanation might be that the mechanism that produces the crystalline 
silicates in the comet forming region is more efficient for small grains than for large 
grains. Both local flash heating events (like shock annealing and lightning) as well 
as radial mixing are more efficient for small grains. For the local production mech- 
anisms this is due to the fact that small grains are more easily heated than large 
grains. In the radial mixing models this is due to the fact that small grains more 
easily couple to the gas, and are thus more easily mixed outwards by, for example, 
turbulent radial mixing. It is, however, unclear if the size dependencies of the var- 
ious models are strong enough to explain the absence of large crystalline silicate 
grains in the coma of Hale-Bopp. 



5 Conclusions 



We have successfully modeled the thermal emission and the degree of linear po- 
larization of radiation scattered by grains in the coma of comet Hale-Bopp. Our 
method has the following important characteristics relative to previous studies. 
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• Both grain size and grain shape effects are taken into account in the calculations 
of the optical properties. 

• The parameters of the best fit model are determined using an objective least 
squares fitting routine. 

• The abundances of the chemical elements observed in interplanetary dust parti- 
cles and meteorites could be used as constraints for the model. 

• The resulting model is consistent with the infrared emission spectrum observed 
in the wavelength range 7-120 jjm and with observations of the degree of lin- 
ear polarization at various phase angles and twelve different wavelengths in the 
optical to near infrared part of the spectrum. 



To model the effects of grain shape on the optical properties, we employed the dis- 
tribution of hollow spheres. In this distribution we average over the volume fraction 
occupied by the central vacuum inclusion while preserving the material volume of 
the particles. We showed that this shape distribution is successful in reproducing 
the observed properties of cometary grains. 



We deduced from the ratios of the strengths of various forsterite features in the 
observed spectrum of Hale-Bopp that the crystalline silicate grains have a volume 
equivalent radius r < 1 jjm. This is much smaller than the typical grain size of the 
other dust components and is in agreement with the sizes of the crystalline silicate 
inclusions found in fluffy interplanetary dust parti cles. The crystalline inclusions in 
these grains are predominantly submicron sized (Bradl ey et al. . 19991) . 



The long wavelength observations showed that most of the mass resides in rela- 
tively large grains. The lack of large crystalline silicate grains in our model thus 
implies that the amount of mass in this component is small. Our best fit model 
has a relative amount of crystalline silicates that is significantly lower than found 
in previous studies of the infrared spectrum. If we consider only the grains with 
a volume equivalent radius smaller than lO/mi, the fraction of the total dust mass 
contained in crystalline silicates is only ~7.5%. The fraction of crystalline silicates 
in the silicate component (both amorphous and crystalline olivine and pyroxene) 
is ~12.5%. This crystallinity can easily be produced by models in which the crys- 
talline silicates are formed close to the Sun by thermal annealing and then mixed 
outwards to the comet forming region (~20-30AU). The crystallinity derived by 
us for comet Hale-Bopp is in agreement with these models assuming an approxi- 
mately nominal model of the pro to solar nebula and a formation of the comet at a 
distance of ~30 AU from the Sun. This crystallinity is also in agreement with that 
found in interplanetary dust particles. 
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